Démonstration des lois de Kepler 
à l’aide du calcul différentiel et intégral 


Introduction 


Johannes Kepler est un astronome allemand qui vécut 
de 1571 à 1630. Né de parents pauvres, il obtint une 
bourse pour étudier le culte luthérien, mais il en fut ex- 
clus un peu plus tard à cause de ses positions religieu- 
ses. Sonrefus d’adhérer au christianisme s’ajouta à ses 
problèmes et après un édit contre les protestants, il se 
réfugia à Prague auprès de l’astronome danois Tycho 
Brahe. Celui-ci le chargea de calculer le mouvement 
de la planète Mars en utilisant les innombrables don- 
nées qu’il avait déjà accumulées. 


Figure 1! 
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Kepler découvrit que l’orbite de Mars était une ellipse 
et publia ses résultats en 1609 dans Astronomia Nova. 
Après plusieurs ennuis dans sa vie privée, où il perdit 
femme et enfants, il retourna en Allemagne et publia 
en 1619 Harmonices Mundi dans lequel il énonce sa 
troisième loi? Voici les fameuses lois de Kepler. 


Première loi : Les planètes décrivent des orbites ellip- 
tiques dont l'un des foyers est le Soleil. 


Deuxième loi : La droite joignant le centre du Soleil 
au centre d'une planète balaie des aires égales en des 
intervalles de temps égaux. 


Troisième loi : Le carré de la période de révolution 
d'une planète est proportionnel au cube du demi 
grand axe de l'ellipse. 


Newion fur le premier à donner une véritable démons- 
tration des lois de Kepler ; cela consacrait en quelque 
sorte ses nouvelles théories sur le calcul différentiel et 
intégral et sur sa loi de la gravitation universelle. 


rielle : la somme des forces appliquées sur un objet de 
masse m est égale à sa masse multipliée par le vecteur 
accélération : 


DÂ-ma. 


Le vecteur accélération se calcule en dérivant deux 
fois le vecteur position : 


A 


Figure 2 


B=x0) F +y] , 


Psx iy], 


ā=x"() i +y Oj. 


1. Démonstration de la seconde loi de Kepler 


Nous allons commencer par démontrer la seconde loi 
de Kepler et ce comme corollaire du théorème sui- 
vant : 


Théorème 1. Si une masse se déplace dans un champ 
de force central, en l'absence de toute autre force, le 
vecteur position décrit des aires égales en des temps 
égaux. 


Démonstration : On place l’origine du système de ré- 
férence au centre de force et on utilise les coordonnées 
polaires. 


Figure 3 


P =r() cosol) T +r singu) j> 
F =[r'Ucos9-rt)sind -00T + [r e)sind + rie) cosd) j> 
ā= [r c0s0-7'sind-0-r'sin8.0 -rcos0-(8) -rsin -0"]-7 
+[r" sind +r'cos0-0'+r'cos0-0'—rsin0-(0"} +rcos0-0']-F 
cos -F +sin0-} et 
le vecteur normal unitaire. 050 -j.On 
peut alors réécrire le vecteur accélération ainsi : 


Le vecteur radial unitaire est 


—sing 


a=[r" 10) ](c0s0-7 +sin0: ÿ+[2r0 +r") (-sin 0-7 +c050-7) 


w a= 


Le-ner Jarro sron: 


Comme la seule force en présence est radiale, la com- 
posante normale est nulle : 


On en déduit que : 


d9 
r? a k où k est une constante. 
it 


Désignons par A l'aire balayée par le vecteur position 
pendant l'intervalle de temps 


SX \ 
SA jo 


Figure 4 
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Un élément d’aire en coordonnées polaires est : 
AA xr: A9-Ar , 


Figure 5 
De telle sorte que : 


A= |frarao 
x 
14, 
A= 5 Je 


Mais comme r2d0 = kdt , il s'ensuit que 


feasta, -1). 


Le vecteur position balaie des aires égales en des 
temps égaux. 


C.Q.F.D. 


Nous voyons que la seule condition utilisée dans cette 
démonstration est que le champ de force soit central et 
qu’il n’y ait pas d’autres forces en présence. 


2. Démonstration de la première loi de Kepler 


Pour démontrer les deux autres lois Kepler nous au- 
rons besoin de la loi de la gravitation universelle de 
Newton : 


Gmm 


F 
P 


En fait, il suffit que la force soit inversement propor- 
tionnelle au carré de la distance. Cette idée était déjà 
dans l’air depuis un certain temps, ce n’est pas Newton 
qui l’a inventée’. Elle s’explique facilement pour un 
champ de force qui rayonne dans l'espace à partir d'un 
point central : comme la surface d’une sphère est 4xr? , 
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la densité des lignes de force est inversement propor- 
tionnelle au carré de la distance. Entre autres, cela est 
vrai pour la gravitation, pour les champs électriques et 
pour l'intensité d'une source lumineuse. 


La première loi de Kepler sera une conséquence du 
théorème suivant : 


Théorème ? Si un objet se déplace dans un champ de 
force central inversement proportionnel au carré de la 
distance entre le centre de l'objet et le centre de force, 
alors la trajectoire suivie par l'objet est une conique. 


Démonstration : Cette démonstration est la suite de la 
démonstration précédente. La farce de gravité est ra- 
diale et dirigée vers le centre ; à partir de la formule du 
vecteur accélération obtenue en (1), nous obtenons : 


dð j ` 
Et, comme r? T , il s’ensuit que : 


CES 


2) 


Comme le but est de trouver la trajectoire décrite par 
l'objet, il faut trouver r en fonction de @ plutôt qu’en 
fonction du temps. De plus, nous effectuerons le chan- 


gement de variable u 


LL 
T il s'ensuit 


En portant cette expression dans l'équation différen- 
tielle (2), nous obtenons : 


d'u K 
6) A 


C’est une équation linéaire d’ordre 2 non-homogène 
de la forme ay" + by'+ cy = c . Ce type d’équation est 


étudié en Calcul III * . On commence par trouver la so- 


du 


lution générale de l'équation homogène W +u=0. 


Pource, on calcule le discriminant du polynôme carac- 
téristique : A =h? — 4ac = —4. Comme sa valeur est 
négative, on aura comme solution 


u= Ke” sin(B0)+ Ke“ cos(0) 


a 


B. 
a’ 
et 


qui donne 0 dans ce cas-ci, 


qui donne 1. On obtient donc 
4(9)= Kı sin0 + K2 cos0 
pour l'équation homogène. Il suffit de trouver une so- 


lution particulière à l'équation non-homogène ct de les 
additionner. Ainsi, 


K 
u(0)= Ki sin0+Kzcos0 +77 


est la solution de l'équation (3). 


Si nous plaçons l'axe des coordonnées polaires de telle 
sorte que r(9) soit minimal lorsque 8 =0, u(8) de- 
vrait être alors maximal, ce qui donne Kı =0 et, puis- 
que w°(0)=-K2 <0, Kı >0 . Ainsi, 


8)= Ki sin0+ Kacos0+ À 


4) = —— #7: 
TES 


Comme nous allons le voir, cela représente l’équation 
d’une conique en coordonnées polaires. 


CQFD. 


3. Équation des coniques en coordonnées 
polaires 

Les coniques peuvent toutes être définies en fonction 

d’un point F (dit le foyer), d’une droite A (dite droite 

directrice) et d’un nombre positif e (appelé excentri- 

cité) : 


dP, F)=e-d(P, A) 


Figure 6 


Dans le cas d’une parabole, e= 1 ; dans le cas d'un” 


ellipse, e= = <1 et dans le cas d’une hyperbole, 


a TE ns P P 
e==> 15. Le cercle estun cas limite où la droite direc- 
a 


trice est à l’infini et l’excentricité est nulle. 
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Trouvons maintenant l'équation d’une conique en 
coordonnées polaires en plaçant le pôle au foyer et en 
prenant l'axe perpendiculaire à la droite directrice. 
d(P,F)=e-d(P, A) 
On voit que 
d(P, F)=r et d(P,A)= p-x= p-rcos0: 
d'où, 


r=e(p-rcos0) = ep—ercosé 


-er 
6 TT Tre-cos0 


L’équation (4) peut se ramener à cette forme : il s’agit 
donc d’une conique. 


© rn0-— + 


— 
# {tt cose) 


F 
k 
-— K 
k? 

1+ g cos? 


2 


Kak 


avec e= . Cela complète la démonstration du 


théorème 2. 


Pour achever la démonstration de la première loi de 
Kepler, examinons les situations possibles. Si un objet 
entre dans un champ gravitationnel quelconque, il 
peut se produire trois cas : si sa vitesse linéaire est suf- 
fisamment grande, il sera dévié par le champ de force 
sur une hyperbole et poursuivra sa route sur une des 
droites asymptotiques. 
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Figure 7 


Par contre, si la vitesse de l’objet est comprise entre 
certaines valeurs, il entrera en orbite sur une ellipse au- 
tour de lastre, sa vilesse déterminera les valeurs de a 
et de b. Dans le cas d’un satellite, cela déterminera 
l’altitude à laquelle il tournera. 


Figure 8 


Finalement, si sa vitesse n’est pas assez grande, il s’é- 
crasera sur l’astre en décrivant une parabole. C’est ce 
qui est arrivé à deux reprises aux satellites de la NASA 
sur Mars 


Figure 9 


Comme la trajectoire des planètes autour du Soleil est 
fermée, on peut conclure qu’il s’agit d’ellipses, ce qui 
termine la démonstration de la première loi de Kepler. 


4. Démonstration de la troisième loi de Kepler 


Démontrons maintenant la troisième loi, en poursui- 
vant la même démarche que dans les deux autres dé- 
monstrations. 


Supposons que la trajectoire soit une ellipse : l’équa- 
tion (6) peut alors s’écrire sous la forme : 


re)= K(L+ ecos) ` 


La distance minimale de 7 sera atteinte lorsque le co- 
sinus atteint son maximum à @ = 0 ; d’où 


k 
Tan = TATA 
K(+e) 

et la valeur maximale sera atteinte lorsque le cosinus 
atteint son minimum à ĝ = x: 


frs 

mn = Kae) 
Soient P le point le plus rapproché du foyer F; et P' 
le plus éloigné ; 


2 F; 
Figure 10 
PF, + P'. 
Tmin + Tmas k? 
a: = s 
2 Ke’) 


Dans une ellipse, l’excentricité est égale à c/a, où c est 
la distance focale et c? = a? - b’? . Ainsi, 


d'où 
k o koa 
b Kb 
P 


Soit T la période d’une révolution complète ; comme 
il a été démontré précédemment, laire balayée par le 
vecteur position entre le temps ñ etle temps #2 est 
donnée par la formule 


k 
A= 34 -h). 
L’aire totale de l’ellipse étant nab , nous obtenons : 


ab=lur, 
2 


K 


ce qui démontre la troisième loi de Kepler. 


On peut faire refaire à des étudiants de Calcul III toutes 
ces démonstrations (en donnant évidemment beau- 
coup d’indications) comme travail de session dans le 
cadre de l’activité d'intégration. C’est un bel exemple 
de lien entre la physique et les mathématiques. En leur 
faisant faire une introduction historique plus détaillée, 
on touchera aussi à l’histoire des sciences. W 
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Notes 


l Photo tirée du site internet http:/www.kepler 
-arc.nasa.gov/johannes.html : elle provient du musée 
autrichien der Sterwarte Kremsmünster. 


? On peut consulter la référence [1] pour en savoir plus 
sur Kepler. 


3 On peut consulter [6], un magnifique livre consacré à 
la vie de Newton. 


* Pour plus de détails sur cette méthode, on peut con- 
sulter la référence [3], pages 307 et suivantes. 


? Le lecteur aura deviné que a représente le demi-axe 
horizontal et b, le demi-axe vertical, dans le cas d’une 
ellipse. Pour les équations des coniques en 
coordonnées cartésiennes et pour une étude plus 
détaillée de celles-ci, on peut consulter la référence 
m. 
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